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ABSTRACT  t 1 
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Consider  independent  random  variables  , X2» • • • ,Xn,  uniformly  f 

\  INSF 

distributed  on  unequal  intervals  (0,a^),  a^  real,  i  =  1,2,...,n.  Let  \. 

S  =  X.)  +  X2  +  • . .  +  Xn.  The  density  of  S  is  shown  to  be  a  linear 
combination  of  ni  B-splines  with  constant  coefficients.  A  probabilistic 
interpretation  of  the  B-spline  as  the  density  of  a  linear  combination  of  order 
statistics  from  the  uniform  distribution  on  (0,1)  is  given.  This 
interpretation  makes  it  possible  to  establish  recurrence  relations  for 
densities  and  moments  of  linear  combinations  of  order  statistics  and  to  give 
asymptotic  results  both  for  B-splines  and  linear  combinations  as  well. 

Examples  and  applications  are  also  discussed. 
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SIGNIFICANCE  AND  EXPLANATION 

B-splines  and  linear  combinations  of  order  statistics  play  an  important 
role  in  the  approximation  of  functions  and  statistics  respectively.  It  is 
useful  to  investigate  the  relation  between  the  two  subjects,  both  from  the 
theoretical  and  computational  point  of  view. 

B-splines  on  equidistant  knots  were  viewed  nondeterministically  by  I .  J . 
Schoenberg  already  in  his  first  work  on  B-splines.  Several  studies  on 
distributions  of  sums  of  independent  uniformly  distributed  random  variables 
have  given  formulae  which  in  fact  are  easily  seen  to  represent  splines. 
However,  a  precise  probabilistic  interpretation  of  the  B-spline  seems  not  to 
be  explicitly  stated  in  the  literature. 

t  r\  f  hr  C  t  y  i  ^ 

!  It  is  shown  'here  that  the  density  of  the  sum  of  n  independent  random 
variables  uniformly  distributed  on  unequal  intervals  is  given  by  a  linear 
combination  of  nl  B-splines  with  constant  coefficients. 

Another  useful  representation  of  the  same  density  is  given  using 
de  Boor's  definition  of  the  discrete  B-spline. 

It  is  also  shown  that  the  B-spline  is  the  density  of  a  linear  combination 
of  order  statistics  from  the  uniform  distribution  on  (0,1).  This  interpre¬ 
tation  of  the  B-spline  allows  one  to  establish  easily  validity  of  recurrence 
relations  for  densities  and  moments  of  linear  combinations  of  order  statistics 
as  well  as  to  state  limit  theorems  both  for  B-splines  and  linear  combinations. 

Examples  given  below  illustrate  how  asymptotic  results  for  B-splines  may 
be  applied  to  linear  combinations  of  order  statistics.  Using  limit  theorems 
of  probability  theory  two  examples  of  Curry  and  Schoenberg  (1966)  for  B- 
splines  are  derived  even  in  somewhat  greater  generality  (see  example  3)1 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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B— SPLINES  AND  LINEAR  COMBINATIONS  OF  UNIFORM  ORDER  STATISTICS 


Z.  G.  Ignatov  and  V.  K.  Kaishev* 


1.  Introduction  and  nanary. 

Let  Xf,X2/>-><Xn  be  independent  random  variables  and  X^  be  uniformly  distributed 

on  the  interval  (O.a^),  aA  >  0,  i  “  1,2,...,n.  Let  S  «  X1  +  X2  +  . ..  +  Xn.  An 

expression  for  the  density  of  S  was  found  by  Olds  (1952),  applying  the  convolution 

formula  for  the  density  of  a  sum  of  two  random  variables  and  an  induction  argument. 

Dempster  and  Kleyle  (1968)  have  given  a  geometric  derivation  of  the  distribution  of  S 

n 

under  the  condition  £  X. /a,  *  1  and  have  pointed  out  the  relation  of  the  latter  to  the 
i-1  1  1 

distribution  of  linear  combinations  of  uniform  order  statistics  on  (0,1).  The  distribution 
of  such  linear  combinations  of  order  statistics  in  more  generality,  covering  the  case  of 
zero  terms  in  the  linear  combination  has  been  obtained  by  Ali  (1969)  and  Weisberg  (1971). 
Cicchitelli  (1976)  approached  the  same  distributional  problem  by  direct  integration.  A 
related  distribution  was  recently  considered  by  Currie  (1981). 

In  the  present  note  the  relation  of  the  above  mentioned  distribution  to  polynomial 
splines  is  emphasized.  Theorem  2.1,  Section  2  gives  an  exact  representation  of  the  density 
of  S  as  a  linear  combination  of  nl  B-splines.  The  B-spline  Mn( x>xQ ,x1 , . . . ,xR) 
introduced  by  Curry  and  Schoenberg  (1966),  of  order  n, 

(degree  n  -  1)  with  knots  xQ  <  x1  <  . . .  <  xn  is  defined  by  the  familiar  divided 


difference  representation 


n  (x  -  x)  + 

<1-1)  M„(X,X0'X1 . Xn)  •  "  ,l0  ,|(7j) -  ' 

where  u>(x)  -  (x  -  xQ)(x  -  x1  )  ...  (x  -  xn>,  (x)+  »  max  {x,o}  . 


•This  work  was  initiated  at  the  Institute  of  Mathematics,  Bulgarian  Academy  of  Sciences  and 
completed  while  the  second  author  was  a  participant  in  the  scientific  exchange  program 
between  the  National  Academy  of  Sciences  of  United  States  and  Bulgarian  Academy  of 
Sciences. 
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Let  zvz2' *  ’  *  ,zn  be  independent  uniformly  distributed  random  variables  on  (0,1). 

Let  *  Z(2)  *  •••  *  Z(n)  be  the  order  statistics  of  Z1 , Z2, • • • . Zn.  It  is  shown  in 

Section  3,  Lemma  3.1  that  the  B-spline  coincides  with  the  density  of  the  linear 
combination  L  =  a^z^j  +  &2Z{2)  +  •••  +  anz(n)*  Recurrence  relations  for  densities  of 
linear  combinations  of  order  statistics  and  their  moments  are  also  given.  In  particular. 
Lemma  3.2  of  Section  3  gives  a  probabilistic  analogue  of  the  well  known  recurrent  relation 
for  the  univariate  B-splines  (c.f.  de  Boor  (1972),  Cox  (1972))  in  terms  of  densities  of 
linear  combinations  of  order  statistics.  On  the  basis  of  Lemma  3.1  and  the  limit  theorems 
6  and  7  Section  5,  of  Curry  and  Schoenberg  (1966)  for  the  B-spline,  two  limit  theorems  for 
the  linear  combination  of  the  order  statistics  Z{1 ) »z(2 ) ' *  *  * ,z(n )  are  stated  (Section  4, 
Theorems  4.1,  4.2).  Theorem  4.3  of  Section  4  gives  a  necessary  and  suficient  condition  for 
convergence  of  the  B-spline  to  1/V2s  exp(-x  /2)  when  n  tends  to  infinity.  Under  this 
condition.  Theorem  4.4  states  an  expansion  for  ( t  jx^  , . . .  )dt .  Both  theorems  4.3 

and  4.4  are  easily  obtained  by  correspondingly  applying  Lemma  3.1  to  the  results  of  Hecker 
(1976)  and  van  zwet  (1979)  for  linear  combinations  of  uniform  order  statistics. 

In  Section  5,  examples  and  applications  of  the  results  of  Section  3  and  4  are  given. 
They  Illustrate  both  the  stochastic  approach  to  studying  asymptotic  behavior  of  fundamental 
spline-functions,  and  the  use  of  spline  theory  techniques  in  probabilistic  problems. 


2.  B- opline 8  and  the  distribution  of  S. 

We  first  state  the  following  Lemma  which  interprets  the  B-spline  stochastically. 

LEMMA  2.1. 

The  conditional  density  function  fg(x/X^/a^  >  X2/a2  >  ...  *  Xn/an>  coincides  Vx 
with  the  B-spline  M^ xjxQ ,xt , . . . ,x„)  having  knots  x0  =  0,  xn  =  a,,  x2  -  a,  +  a2,..., 

xn=a1  +a2+  +an* 

PROOF. 

As  assumed  above  the  random  vector  is  uniformly  distributed  in  the  n- 

dimensional  hypercube  (0,a^)  x  (0,82)  x  •••  x  (0#*n)  3  A*  The  random  event 

|x  /a  >  X  /a  >  . • .  >  X  /a  }  *  0  c  ft  is  an  n-simplex  denoted  by  a  with  vertices 
1  1  1  2  2  n  n 1 

(0,0, ...,0),  (a-j,0,0,  •••#  0),  ( a  ^  ,a2  #  0  #  0  #  •  •  •  #  0  )  #  •••#  (*i  •  “  *  *  #an^  *  Then  the 

probability  of  0  is 

n 

-  rr«i 

,  .  vol (0)  n  i-1  1 

p*°>  vol(fi)  n  nl 

TT  *i 

i-1 

Now,  by  definition  fg(x/X1/a1  >  X2/a2  >  ...  > 


lim  - 
h*0 

nl 


?{ (X1  +  Xj  +  . . .  +  Xn  e  (x,x  +  h))  n  (Xj/a,  *  x2/«2  *  >  ^^n^ 

h  Plx1/a1  >  X2/a2  >  ...  >  Xn/an] 


TT'* 

1-1 


lim 

h*0 


vol  (a(x,x  +h)  n  c) 
h 


where  a(x,  x  +  h)  =  ( (u1  , . . .  ,un  )  :  +  u2  +  .  •  .  +  un  e  [x,x  +  h]  }  . 

If  the  hyperplane  u,  +  u2  +  . . •  +  u„  -  x  is  denoted  by  a(x)  then 


lim 

h*0 


vol(a<x,  x  +  h)  n.oj  _  _L  vol  (a(x)  n  a)  , 

h  rr  n_1 


where  voln_t  is  the  (n  -  1)  -  dimensional  volume. 
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Thus  fgtx/X/s,  >  Xj/Sj  >  >  *n/an) 


,  vol  ,  (a(x)  no) 
m  n- 1 


TT  *i 

i-i 


Let  us  assume  that  vol(fi)  -  |  [  a.  «  ni.  Then 

i-1 


(2.1)  f.(xA./a  >  X  /a  >  ...  >  X  /a  )  -  —  vol  .(o(x)  n  o) 

S  I'  2  2  n  n  ^  n- 1 


and  vol(o)  -  1.  Now  we  can  apply  Theorem  2  of  Curry  and  Schoenberg  (1966)  which  states: 
The  fundaaental  spline  function  «„( xi x0, x, , . . .  ,xn )  is  the  linear  density  function 
obtained  by  projecting  orthogonally  onto  the  x-axis  the  volume  of  an  n-dimensional  simplex 
an,  of  volume  unity,  so  located  that  lta  n  ♦  1  vertices  project  orthogonally  into  the 
points  x0,x.|, . . .  ,xn  of  the  x-axia,  respectively. 

Let  us  project  orthogonally  the  aiaplex  o  onto  the  g-axis  with  a  parametric  equation 
u1  «  t/Zn  ,  Uj  -  t/Zn  ,  ....  un  -  t/Zn  ,  t  c  R.  If  the  coordinate  system  0,g  is 
introduced  onto  the  g-axis  so  that  the  point  (t/Zn  ,  t/Zn,  ...,  t//n) 6  g  has  the 
coordinate  t  then  the  orthogonal  projection  of  the  vertex  A1(a1 ,a2 . . . ,ai,0,0, . . . ,0 )  of 
the  simplex  o  onto  the  g-axis  will  be  the  point  with  the  coordinate 

(ai  +  a2  ♦  ...  +  a^/Zn  ,  i  -  1,2,  ...,n,  with  respect  to  0,g.  Then  according  to  the 
above  theorem  of  Curry  and  Schoenberg 

vol^^atx)  n  a)  -  Mn(x»0,  a^/Zn,  («1  +  a2  )//n, . . . ,  (a1  +  a2...  +  an)/Zn)  . 


Now  from  (1.1),  setting  a„  =  0  we  have  vol  , (a(x)  n  o) 

u  n-i 

n  ((aQ  +  8^  ...  +  a  .)/Zn  -  x/Zn)”  1 

*  n  l  - —  Z - — 

j”0  -pj  <(aQ  +  a1  +  ...  +  a..)/Zn  -  (aQ  +  a1  +  ...  +  a^/Zn) 

i-0 

-  Zn  VxiO.a,,*,  +  a2 '  a1  +  a2  +  +  an>  ' 

n 

From  (2.1)  for  the  case  |  |  a^  »  nl  we  have 

i-1 

fs(xAi/a1  >  ...  >  Xn/en )  —  Mn(x/0,a1  ,a1  +  a 2,  •  ••,  a^  +  a2  +  ...  +  8^ ) • 


-4- 


tim 


n 

Let  vol(Q)  -  |  |  a,  *  k  *  nl,  k  >  0.  By  the  simple  transformation 

i-1  x 

Xi  "  d  —  ^nl/k  ,  i  «  1,2, ...,n  we  obtain  the  set  of  random  variables  X.,  X_,...,X 

1  2  n 

uniformly  distributed  on  (0,^),  bt  -  ajd,  i  -  1,2,...,n  and  such  that 
n 

vol(fl)  «  "J  [  b.  -  nl. 
i-1 

Then 

(2.2)  f  (x/X/b.  >  X/b.  >  ...  >  X/b  )  - 

S  1  2  n 

Mn(x»0,b1 ,b^  +  b2,-..,b1  +  b2  +  ...  +  bn)  . 

Changing  back  to  the  variables  X1,X2,...<Xn  on  the  left  side  of  (2.2)  and  applying  (1.1) 
to  the  right  side  we  have 

fg'x/X.,/*,  >  X2/a2  >  ...  >  Xn/an)  - 

Mb'*'0*4,.*)  ♦  +  *2  +  **•  +  V' 

where  x  »  ^  ,  proving  the  assertion  of  Lemma  2.1. 

Remark:  For  clarity  and  simpler  notation  we  have  confined  ourselves  to  the  case  a^  >0, 
i  -  1|...,n.  Since  the  divided  difference  (1.1)  is  a  symmetric  function  of  its  arguments 
ie.,  it  is  independent  on  the  order  of  the  points  xQ,x1,...,xn  in  the  argument  list  it 
can  be  easily  seen  that  the  above  Lemma  is  true  for  a^  real. 

We  are  now  ready  to  prove 
THEOREM  2.1. 

Let  Xj ,X2, . . . ,Xn  be  Independent  random  variables  and  X^  be  uniformly  distributed 
on  the  interval  (O.a^,  aA  -  real,  i  -  1,2, ...,n.  Let  S  -  X,  +  X2  +  ...  +  xn. 

Then 


(2.5,  fs<x)-^ 


nl  l  M  (xjOja^  ,a  +  a.  ,...,a.  +  a.  +  • • .  +  a.  )  , 

‘V1! . V^n  1  1  2  1  2 


where  Qn  are  the  permutations  of  (1,2,...,n). 
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PROOF. 


Vx) 


P  (X.  /a.  >  X.  /a.  >  ...  >  X  /a.  } 

1  l2  x2  n  n 


. vc«n 


fs(x/Xi  /a±  >  XL  /aA  >  ...  >  XA  /at  ) 

112  2  n  n 


1  -f  f,(xA,  /a.  >  X  /a  >  ...  >  X  /a  ) 

<V*a . V  =  «„  '  ’  2  2 


1 

nl 


1  M  (xiO,a  ,a  +  a  , 

. in>cQn  1  1  2 


'V  +  *1,  +  +  “i  ) 

1  2  n 


by  Lemma  2.1. 

Another  interesting  representation  of  f  _  (x)  arises  if  the  notion  of  discrete  B- 

splines  as  viewed  by  C.  de  Boor  (1976)  is  assumed. 

Let  0  <  x1  <  . . .  <  Xjj,  aj*xj,...,an-xn-  x^j  and  ^(i^  , . . .  ,in)  • 

»  {o,a.  #a.  +  a.  ,...,a.  +  ...  +  a.  }  where  (i1>...,in)  c  Qn.  We  shall  also  use  the 
1  1  2  1  n 

notation  T(i  , . . . ,i  )  «  { , .  ....... t  .  .}  in  which  for  the  sake  of 

1  n  0f(i  /•••,!  )  n  /  ( i 4 1  • • • J 

in  in 

brevity  we  will  further  drop  the  second  subscript  and  write  {tg,...,Tn}. 

t 

P 

ordering  the  set 


Let  t^  ■  (t  ,t  ,  ...,t  },  p  “  \  ("J  be  the  strictly  increasing  sequence  obtained  by 

p  i-0 


{o,a.  ,a.  +  a.  ,...,a.  +  a.  +  ...  +  a.  ,...,a.  +  a.  +  . . .  +  a.  }  , 

1  *1  X2  A1  X2  Xj  l1  X2  n 

where  (ij,...,ij)  c  c”,  are  the  combinations  of  j  from  n,  j  «  1,...n. 

The  B-splines  of  order  n  for  the  knot  sequence  are  given  by 

Mi,n,t<x)  “  n^ti'",,ti+n^'  “  x)+_1'  1  "1****»« * '  where  [t± , . . .  ,ti+J  denotes  a 

n-th  order  divided  difference  and  q  “  p-n. 

Since  H_(xiO,  a.  ,a.  +  a.  +  a.  +  . . .  +  a .  )  of  (2.3)  is  also  a  spline 

"  1 1  l2  x2  *1  x2 

with  knots  _t.  it  can  be  represented  as  a  linear  combination  of  the  B-splines  Mi,n,t^x^' 
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1 


i  «  1,...q  by  the  Curry-Schoenberg  (1966)  theorem  (see  p.  113).  The  non-negative 
coefficients  of  such  a  linear  combination  termed  by  C.  de  Boor  (1976)  discrete  B-splines 


are  defined  as 


(2.4)  8 


,n,T(i1 . i„),t(i)  =  <ti+n-ti)tTj . Tj+nll‘  -  fci+lV  '•  -  ti+„-i>+  * 


Remark  1:  When  j  =  0  we  simply  write  8 


Remark  2 :  When  n  =  1 ,  (2.4)  reads 


j  /  1  r 7 (i j  /  •  •  •  / i^  ) / 1 


n  rlS 


iJ  -  <ti+rti>[Tj'Tj+1]('  -  V* 


VO  ,  x  <  t. 


Thus  the  density  f s ( x )  is  expressed  in  terms  of  discrete  B-splines  as 

.  q 

(2.5)  f  (X)  =  l  l  S  (i)M.  (x)  . 

1-1  (i  ,...,i  )cq  n,l<ir-‘*,in)'i  x'n'i. 

1  n  n 

Note  that  if  tj  <  x  <  fcj+i  '3  =  n,...,q-1  then 

fS(x)=^7  S  i  0„  T(i  i  )  t'A)Mi  n  t,X)  ' 

i-j-n+1  (i  ....,i  )c0  n'I<V"V'-  X'n'£- 

1  n  n 

since  other  vanish  identically.  This  means  that  the  evaluation  of  fg(x) 

for  any  x  e  (xo/*nl  requires  summation  of  no  more  than  n  of  the  q  terms  on  the 
righthand  side  of  (2.5). 

The  discrete  B-splines  in  (2.5)  can  be  determined  from  the  recurrence  expression 

®  j  »n,T(  i  , . . .  ,i  ),t*x*  aj,i,n^T  j-t-n  fci+n-1  j  + 1  ,n-1  * x  ' 

—  in  — 

(2.6)  +  (ti*n-1  ’  Tj)ej,n-1(i)l  •  where 

“j.i.n  “  (titn  -  V/l<Vn  '  V'Wl  *  V1  * 


Remark:  For  a  proof  of  (2.6)  with  differently  normalized  discrete  B-splines  see 

Jia  (1983). 


»  •  v  *  »  *  *.  m'  «"  ,»  *  •  •  ’  .  -  .  ■ 

jt  ’  ■  *  ,  •  .  *  .  »  .  i  S  *.  V  •  >  •  >  *  • 

/V'  V v  V  .v  ■.*  . 

““  -  ■  • - •-  •  ■  *  fc  1 m  »  «  .  »  ,  »   s  ,*  » 
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3.  B- splines  and  order  statistics. 

We  first  state  the  following  Lemma,  providing  a  probabilistic  interpretation  of  the  B- 
spline  in  terms  of  order  statistics.  It  should  be  noted  that  this  result  easily  follows 
differentiating  the  distribution  formula  given  by  Dempster  and  Kleyle  (1968)  and  Ali 
(1969).  For  the  sake  of  consistency  with  the  rest  of  this  note  another  proof  is  given 
below. 


LEMMA  3.1. 

Let  Z1,Z2,...,Zn 
variables  and  let  z(i) 
be  fixed  real  numbers. 

a2Z(n)  +*"+  anz(n)  is 


be  n  independent  uniformly  distributed  on  (0,1)  random 
>  Zj2j  >...>  2 j n j  be  their  order  statistics.  Let  a.j,a2,. 
Then  the  density  of  the  linear  combination  L  *  a^Z^j  + 
fL<x)  =  Mj^xjO,  a.|,a.,+a2,  ••  .,a.,+a2  +  .  •  .+an). 


PROOF. 

From  Lemma  2.1  we  have 


I  7  +.  *7  +  x.  .  (X  /  Z.  >  Z.  ^*  .  •  ^  Z  ) 
i1Z1+a2Z2+.. .+anZn  1  2  n 

*  Mn(xj0,a1/a1+a2, . . .  ,a.|+a2+. .  .+an) 


It  can  be  easily  seen  that 


(3.1) 

where 


^a  Z  +a  Z  +. .  .+a  Z  <x  '  Z1  >  Z2  Zn  > 

112  2  n  n 


=  lim 
h+0 


P{ (a.Z,ta0Z_+. . .+a  Z  e  (x,x+h)>  n  (z,  >  Z-  >...>  Z  )} 
IT  z  2  nn  i  l  n 


h  P{zt  >  Z2  >...>  Zn) 


lim  —  f  ...  J  1  du,...du  *  lim  /  ...  f  n!  du, ...du  , 

^  h  1  n  v._.n  h  ;  ;  1  n 


b+0 


h+0 


A,nA2 


4  *  { (u^  ,Uj  , . . .  ,un  )  :  1>u1>u2>...>un>0  ,  x<a1u1+a2u2  +  .  .  .+anun<x+h}  , 
A,  =  { (u1  ,u2 , .  . .  ,un  )  :  1>ut>.. .>un>0}  , 

A2  =  {  (u1  ,u2  , . .  .  ,un  )  :  x<a1u1+a2u2  +  .  .  .+anun  <  x+h) 
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II 


l 


•/%’  %  _ 


.  \ 


Denote 


i  ,  (u1 ,u2 , . . . ,u  )  e  i1 


i  v°  «  <Ui,u2,...,un)  e 


Then  (3.1)  yields 


lim  1  /  ...  /  nl  du, . . .dun  . 


h-M)  A2  1 

Since  the  density  of  the  random  vector  ^z( 1 ) >z( 2 ) ' *  *  * ,z(n) ^  is  *"A  n!  (c,f’ 


S.  Karlin  and  H.  Taylor,  (1981),  p.  102) 


(3.3)  £  J  ...  /  CA  nl  dUl...dun 


PU1Z(1)+a2Z(2)  +  -”+anZ(n)  e  (X'X+h)} 


From  (3.3)  replacing  in  (3.2)  we  finally  obtain 


lim /  ...  /  5,  nl  du,du,  ...du 
.  ^  h  1  .  1  d.  12  n 


h+0  d„  1 

2 

pU1Z(1)+a2Z(2)+"-+anZ(n)  e  <X'X+h))  .  .  . 

*  lim - - -  “  fT(x) 

,  *  h  i. 

h+0 


Since  the  vector  {Z( , } ,Z(2 j , . . . »2<n) '  haa  a  density»  the  random  variable 

a^j  1  )+»2z[  2  )*’  *  *+anz(n)  a^so  haa  a  density,  denoted  by  f^fx)  in  (3.4).  That  completes 

the  proof  of  Lemma  3.1. 

Remark  IS  It  can  be  easily  seen  by  a  linear  transformation  that  if  xQ  <  x,  <...<  xn,  a1  = 

x.,-x0,  a2  ”  Xj-x,,...,  an  -  Xjj-Xj,..,  and  Z1,Z2,...,Zn  are  independent  random  variables 

uniformly  distributed  in  the  interval  (Xg/(xn-Xg)  ,  xQ/(xn-x0)  +  1)  then  the  linear 

combination  L  -  a.Z, ,+a_Z,_+.  .  .+a  Z  has  the  density  f  (x)  =  M^(x;xn,x. , . . . ,x  ) . 

l(l)2(2)  n  (n)  ~  11  ui  u 

Since  «  Zi  +  x0/(xn-Xg)  ,  i  “  1,2,  ...,n  and  Z(ij  *  z(i)  +  *0/,*xn~x0*'  1  =  1,2, ...,n, 

f~(x)  »  fL+x  (x)  -  (^(xixg.x,, . . . ,xn)  holds. 

L  0 

Remark  2.  If  some  a^,  i  ”  1,...,n  are  zero  i.e.  the  B-spline  has  multiple  knots  than  the 
above  Lemma  still  holds.  The  righthand  side  of  formula  (1.1)  must  be  replaced  by  the 
expression  for  divided  differences  for  multiple  arguments  (c.f.  Shumaker  (1981)  p.  119). 

For  stable  computation  of  fL(x)  however  the  recurrence  relations  given  next  should  be 
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Application 


The  circular  serial  correlation  coefficient  with  lag  A  is  given  by 

N  _  _  N  _  2 

(5.1)  aB„  -  (Xj-X)  (Xj+A~X)/  ^  (X^X) 

where  Xg+j  -  Xj  and  A  <  sample  size  N. 

It  may  be  noted  that  the  marginal  density  of  in  the  case  Xj,  1 , . . .  ,N  i.i.d. 

N(0,1)  and  N-odd  is  given  by  the  B-spline  Mn* x> ac1 »ac2' * • * 'Acn+1 ^ '  where 
AC1  >  aC2  >•••>  ACn+1  are  the  distinct  latent  roots  of  the  matrix  of  the  quadratic  form  in 
the  nuswrator  of  (5.1),  n+1  “  (N-1)/2.  This  follows  differentiating  and  slightly 
rewriting  formula  (1.1)  of  Dempster  and  Kleyle  (1968)  which  gives  the  marginal  distribution 

of  aV 
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L®t  T2n  “  y'3y/n  <zi,2n  +*”  +  Z2n,2n)  “  ^3n*  Then'  for  the  density  of  T2n  we  have 


,,  X,  +. . .+  X„ 

h”'*  --in'  2" 


<  x)  =  f  (x) 
2n 


*  M2n(x,-n  ,,3/n'  (-n+1>  ^Vn,  (-n+2)  /3/n  .....  n  /3/n) 

Then  according  to  the  LCLT 

li«  sup  M  <x,-n/37K,  (-n+1  )/57n,  <-n+2)/37^  .....  n/37£)  -  —  exp(-x2) 

n-M.  -~<x<»  1  " 

EXAMPLE  5. 

T2k-1  “  _  T“  +  2k  7T3  Zl,2k-1  +  2k  3^5  Z2,2k-1  +  2k  5^7  Z3,2k-1 

+  ‘"  +  2k  ( 2k-3 )  ( 2k- 1 )  Zk-1,2k-1  +  2k  (  2k- 1 )  Zk,2k-1  + 

2k  (2k-3)(2k-1)  Zk+1,2k-1  +  2k  (2k-5) (2k-3)  Zk+2,2k-1  +*‘*+  2k  TTJ  Z2k-1,2k-1 


Then ,  by  Lemma  3.1 


^T2k«i<x)  ”  Mn(X,~  '  -f 


2k  2k  2k  2k . 

2k- 1  '  2k- 1  '  2k^3  '*•" 


and  by  Example  5  of  Curry  and  Schoenberg  (1966) 


lim  sup  f  (x) - !— — 

k—  — <x<-  I  T2k-1  w  cosh  x 


EXAMPLE  6. 


Let  T 


2k- 1 


log  k  +  +  2  log  k)  ZL  +  (^L_  -  £*)z 

k  k,2k-1  'k-1  k  '*k+1,2k-1 


2k  2k. 


+  (k-2  "  k-1,Zk+2.2k-1  +,*’+  "  T~)Z2 k-1,2k-1  ”  ' 2  lo9  k  +  2(1  +  log  k)Zk  + 


♦  2k 


1 


k ( k- 1 )  k+1,2k-1 

the  density 


+  2k 


1 


(k-1)(k-2)  k+2,2k-1 


+. ..+  2k 


■jjTY  Z2k-1,2k-1*  Then  by  Lemma  3.1 


fT2k_1(x)  “  M2k-i(x,“2  109  k'-*"  ~2  109  k'  ir  '  fh-  T~ 


)  and 


lim  sup  |  f 
k-x»  -*<x<» 


(x>  -  exp(-exp(-x)  -  x) 

2k- 1 


by  Example  6  of  Curry  and  Schoenberg  (1966). 
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f  exp(-x)  ,  x  >  0 

f_  (x)  converges  uniformly  to  X(x)  -  < 

n  v.  0  ,  x  <  0 

for  x  outside  an  arbitrarily  small  neighborhood  of  the  point  x  *  0. 

EXAMPLE  3. 

Let  Tjjj.)  *  -/k  +  2/k  2k-1‘  ®*en»  be®“«  3.1 

— X 

f  _  ( x )  *  H*.  .  ( x  j  -/k  i  *  •  ■  .  -/k  ,  /k,  ..  .  ,/k) 

*2k-1 

converges  uniformly  for  all  real  x  to  exp (-x2)//*. 

More  generally,  by  a  theorem  of  Mosteller  (c.f.  David  (1981)  p.  255)  it  is  known  that 
if  0  <  p  <  1  and 

T  -  -/SE  +  - 5 -  2 

"  ^  /np~('i-p)  n‘tnl>J'n 

then  TR  converges  weakly  to  the  random  variable  T  normally  distributed  with  ET  *  0 
and  DT  -  1.  Hence,  by  Lenma  3.1 

n-  [np]  [np]+1 

/*_  M(u,-/HTT^.,  -nar£L~7^^i1ZE)_)du  _L  f*_  exp(-u2/2)du 

-  n  /1-p  ^1-p  /np(i-p)  n*~/2W 


n  /npTT^pT  /np(  1-p)  n+**  /2tT 

for  every  x  e  R1 . 

EXAMPLE  4. 

According  to  the  local  central  limit  theorem  (LCLT)  if  Xn,  n  -  1,2,...  are 

2 

independent  with  a  common  distribution  F  having  mean  0  and  variance  0  and  if  F  has 
a  bounded  density,  then 


lim  sup 

4-  P(/n  X  <  x) 

dx  n 

- —  exp(-x2/202)  1  - 

n-*»  -°'<x<“  ' 

/2»  ' 

(Zn  -  1/2). 

Then 

j 

f  ’ 

x  >  1/2/6 

F(x)  -  t 

x//6  ♦  1/2  , 

-1/2/6  <  x  <  1/2/6 

t  o 

x  <  -1/2/6 

and  the  conditions  of  the  LCLT  hold  with  a  m  1/2.  Thus  we  have 

X1+...+X2n  _  _ 

/2^  _  -  /V^  <Z12n  +...+  Z2n  2n)  -  /3n  . 


/V. 


and  applications. 


The  following  examples  illustrate  how,  on  the  basis  of  Lemma  3.1,  providing  a 


probabilistic  interpretation  of  the  B-spline,  asymptotic  results  for  B-splines  (Curry  and 


Schoenberg  (1966),  Section  9)  may  be  applied  to  linear  combinations  of  order  statistics. 


It  is  also  shown  (examples  3,  4)  that  using  limit  theorems  of  probability  theory  the 


results  of  Curry  and  Schoenberg  (1966)  (examples  3,  4)  can  be  derived,  even  in  slightly 


more  generality  (example  3). 


As  above,  here  Z,  „  >  Z,  _>...>  Z_  „  is  the  order  statistics  from  the  uniform 
i,n  2,n  n,n 


distribution  on  (0,1). 


EXAMPLE  1. 


Let  Xq  <  x1  <  . . .  be  a  convergent  sequence  of  reals  with  lim  xR  *  B.  Then  the 

n-H» 


sequence  of  random  variables  T„  -  xfl  +  (x^XjlZ,^  +  <x2-x1)z2,n  +‘**+  (xn~xn-1 )zn,n 


converges  weakly  to  the  constant  B. 


According  to  Lemma  3.1  and  Lemma  6  of  Curry  and  Schoenberg  (1966)  for  the  density 


fT  (x)  of  the  random  variable  T_  we  have 
n 


/!.  fT  <x)(1  +  xs/(n+1))"n-1dx  -  Mn(x)(1  +  xs/ ( n+ 1 ) )  ' 


TT  <1  +  x.s/(n+1)) 

i*=0 


It  can  be  easily  seen  that 


TT<1  +  x  s/(n+1)) 
i-0 


Hence,  by  Lemma  4  of  Curry  and  Schoenberg  (1966)  TR  converges  weakly  to  B  and 


therefore  Tn  converges  to  B  in  probability. 


EXAMPLE  2. 


Let  Tn  -  nZn  n.  Then,  by  Lemma  3.1  the  density  of  the  random  variable  Tn  fTn*x* 
Mn(x;0,0, . . . ,0,n)  and  hence,  by  Theorem  7  of  Curry  and  Schoenberg 


w, 


I 


I 

I 

I 

I 

The  remainder  ten  in  the  lest  expansion  is  given  by  van  Zwet  (1979)  to  be  typically 
of  order  n-3/2.  j 


Me  finally  remark  that  the  Berry-Essenn  bound  of  Corollary  1  can  be  sharpened  by  means 
of  Corollary  2  to  obtain 


•  v.  •?.  o-.  t.  r:  rr  y.; 


T»  1  ».'■ 


Following  van  Zwet  (1979)  we  develop  the  expansion  of  Theorem  4.4  for  the  special 
cases  m  »  3  and  m  -  5.  For  m  «  3  it  reduces  to  a  bound  of  Berry-Essenn  type, 
typically  of  order  n~  1/2  . 

COROLLARY  1. 

There  exists  a  constant  C  such  that  for  every  n  -  1,2,...  and  every 
n  2 

aj(n) , . . .  ,an(n)  with  'i  a.(n)  t1  0, 

j-1  3 

|  Mn(t»y0.n'yl,n . *n,n)dt  "  *(x)|  <  C  I  l*a  ^TnVl)  (n+2)  |  3  . 

X  I  jxQ  J 

now  a  *  5.  Define 


i  n  0 

Fn(x)  -  *<x)  -  ♦(x)l-  l  ly3  // (n+1) (n+2) ) 3  .  H  (x)  + 
j-0  3'n  2 

n 


+  ty j ^ „/^ ( n+ 1 ) < n+2 ) ] 4  -  -1-}  H3<x)  + 


2 

+  la  {  L  ty  .  /<t <n+T)  (n*2)  J 3  j  H,(X)1  . 

j-0  3'  5 


COROLLARY  2. 


There  exists  a  constant  C  such  that  for  every  n  -  1,2,...  and  every  a, (n) , . . . ,an(n) 
n  2 

with  l  a  (n)  ?  0, 
j-1  3 


8"P  I  Mn(t'*0.n'*1,n . *n,n)dt  " 

n  _ 

-  Fn(x)  j  4  C  IJyj>n/V(n-M)(n+2)| 


Corollaries  1  and  2  follow  by  Lemma  3.1  and  Corolaries  1  and  2  of  van  Zwet  (1979). 
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v,  _< x )  -  X  -  X  -  o  x 

j,n  3,n  n  n 


n  ,  -  Va 

W  (x)  -  V  ( x ) |  l  ( x ) } 

J»n  j»n 


f  -2-,V/2  (y,  -  x){(n+1)(n+2)}‘ 1/2  , 

l  n+2+x2J  3'n 


5(X)  -  -  l  W^U)  "  ~~~  1  2^  ^ 


j-0 


.  n+2+x 


It  i.8  easy  to  check  that  £  wj,n*x^  *  1*  For  integer  a  >  3  and  real  x  and  z,  let 


a-1 


Gm#n(x,x,  -  «(x)  -  ♦<*>!*  h(E(cv1)(«,  •  rt^-  «5fn(x,} k  , 


where  1  denotes  summation  over  all  nonnegative  integers  v 


'Vi  with 


1  <  l  (k-2)v  <  m-3  . 

k»3  * 


THEOREM  4.4. 


For  every  integer  u  >  3  there  exists  a  constant  such  that  for  every  n  »  1,2,.. 

2, 


and  every  a^Cn) , . . .  ,an(n)  with  t  aj(n)  +  0, 

sup  j  /!-  Mn(t,y0fn,y1n,...,yn>n)dt 


Gm.n<*<*,'*>  1  <  C»  iQ  |yj>n//(n+D(n+2,|  . 


Follows  by  Lemma  3.1  and  the  theorem  of  van  Zwet  (1979)  establishing  the  Edgeworth 
expansion  for  the  distribution  function  of  Tn. 

Since  we  require  the  condition  (4.4)  to  be  fulfilled  the  remainder  in  the  above 
expansion  tends  to  zero  for  every  m  >  3  and  for  increasing  m.  The  rates  will  depend  of 


course  on  the  triangular  array  (a^(n)}  related  to  the  knots  of  the  B-spline.  The 
remainder  approximately  behave  as  n  ^2  (m-2)  (Bee  van  zwet  (1979)). 
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J 

* 


*3 

•J 


W 

l 


S 


V  V  T"  :7\r 


<  ■  »  i  ■  i  *  v 1  .  »  .  *  .  '■  » v  v  ^  v'  ;  *  - 


yJ,n  -  (xj,n  "  xn,/0n'  3  “  0 . . 


Since  E  *j(n)  ¥  0  iaplie*  that  ¥  0  the  yj#n* 


s  are  well  defined  with 


rj.n 


0  and  £ 


yj,n 


(n+1 )  (n+2) . 


THEOREM  4.3. 


li»  »up  M  (x>y  ,y.  ,...,y  ) - 1-  exp(-x2/2)  -  0 

n+-  — <x<-  n  °'n  1'“  n'n  m 


if  and  only  if 
(4.4) 


lin  [  max  |  y.  //(n+1 )  ( n+2 )  |  ]  »  0 

n+-  0<  j<n  ;,'n 


Applying  Lenina  3.1  to  the  theorem  of  Hecker  (1976)  which  gives  a  necessary  and 
sufficient  condition  for  the  asymptotic  normality  of  the  linear  combination  Tn  -  y0>n  + 

<yi,n-yo,n>21,n  +  <*2,^*1, n>z2,n  +  -”  +  ‘yn.n-yn-l.n^.n  we  e*8ily  obtain  th«  «»»«rtion 
of  Theorem  4.3. 

Under  the  condition  (4.4)  we  shall  derive  an  expansion  for  the  spline  distribution 
function 

/l.  «  (tiy0  ..y,  „""'yn  >dt 

n  u,n  i,n  nfn 

using  a  theorem  of  van  Zwet  (1979). 

Define 

*<x>  -  ♦(t)dt  , 

♦  (x)  -  exp(-x2/2), 

/2’  r 

Hj.  -  the  Hermits  polynomial  of  degree  r  i.e.  — — — —  “  {-1)r$(x)H  (x).  For  j  “  0,...,n 

dxr 

and  real  x  let 
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/  f 


THEOREM  4.2 


Let  the  sequence  of  random  variables  Tn  ■  a^(n)Z^  )  +  a2(n,z(2)  +**,+ 
an(n)Zjn j+an+1(n)  converge  In  distribution  to  a  random  variable  with  a  distribution 


function  F(x),  and  let 


tm  -sx  .  ,  a 

J_ o,  e  dF(x)  ?  e 


Hence  dF(x)  -  A(x)dx,  where  A(x)  is  a  Polya  frequency  function  and 


lim  fT  (x)  -  A(x) 
n+**  n 


uniformly  for  all  x  provided  that  A(x)  is  not  of  the  form  (4.2).  Here  f_  (x)  denotes 


the  density  function  of  the  random  variable  Tn« 


For  A(x)  «  A.(x),  again  (4.3)  holds  uniformly  for  x  outside  an  arbitrarily  small 


neighborhood  of  the  point  x  “  6-6^. 


If  4(a)  -  exp(6s)  as  in  case  i)  then 


lim  f  (x)  -  0 
n—  " 


uniformly  in  x  outside  an  arbitrarily  small  neighborhood  of  the  point  x  ■  6 . 


The  proofs  of  Theorems  4.1  and  4.2  directly  follow  from  theorems  6  and  7  of  Curry  and 


Schoenberg  (1966)  respectively,  applying  Lemma  3.1. 


He  shall  now  use  two  theorems  of  asymptotic  statistics  to  establish  easily  results  for 


B-splines . 


Let  a1(n),...,an(n)  be  a  real  sequence  with  x0  n  "  0,  x1n  -  8,(0),  x2>n  “  a^n) 


11  2 

(n) , .  • .  .Xjj -  a,(n)  +...+  a„(n)  and  I  a^(n)  +  0  for  every  n. 


Define 


1  P 

Kn  "  n+1  xj,n  ' 


l  l*.  -  *  >2 

2  j-0  j,n 


( n+ 1 ) <  n+2 ) 


L-‘.l 


4.  Limit  theori 


Let  us  first  state  two  limit  theorems  concerning  linear  combinations  of  order 
statistics.  For  the  purpose  the  class  of  Polya  distribution  functions  F(x)  will  be 


required,  these  are  the  distribution  functions  having  a  bilateral  Laplace  transform  of  the 


(4.1)  e  8  dF( x)  -  1/*(s)  , 

« 

2  — ■  o 

where  *(s)  »  exp(-ys  +  6s)  [ [  ( 1+6 . s)exp( -6  s)  ,  Y>0,  6,6  real  ,  E  6  <  «•  (see  Curry 

i-1  1  1  1  1 

and  Schoenberg  (1966),  p.  91). 

THEOREM  4.1. 


If  {a1(n),a2(n), . ..,an(n),an+1(n)}  is  a  sequence  of  constants  such  that  the  sequence 
of  linear  combinations  of  order  statistics  a.,  (n)Z^  j+a2(n)Z(  2  j  +  . .  .+an(n)Z^n  j+an+1  (n) 
converges  weakly  to  a  random  variable  W  with  a  distribution  function  F(x),  then  F(x) 


is  a  Polya  distribution  function.  Conversely  if  W  is  a  random  variable  with  a  Polya 
distribution  function  F(x)  then  there  exists  an  appropriate  sequence  of  constants 
a1(n),a2(n),...,an(n),an+1(n)  such  that  a1  (n)Z(  ,  j+a2(n)Z(2  )+•  ..ta^n  )Z(n)+an+1  (n) 
converges  weakly  to  W. 

2 

It  is  known  that  exp ( fa  )  is  the  bilateral  Laplace  transform  of  the  random 

variable  W  normally  distributed  with  EW  «  0  and  D(W)  •  2y.  Consequently  *(s)  1 

(c.f.  (4.1))  is  the  bilateral  Laplace  transform  of  the  distribution  of  the  sum  of 

independent  random  variables  U  -  6  +  /2yw  +  6.^  -  61  +  62V2  -  62+. ..,  where  is 

exponentially  distributed  with  EV^  •  1,  i  «  1,2,...  ,  W  ~  N{0,1),  and  6  real  constant. 

i)  If  0  ■  Y  “  61  -  62  “.. .  then  the  bilateral  Laplace  transform  of  the  random 

variable  U  is  *(s)  -  exp(6s)  and  no  density  function  exists. 

ii)  If  Y  «  0  and  only  one  of  the  constants  6A  is  nonzero,  i  »  1,2,...  then  the 

random  variable  U  has  a  density  function 

{•j— exp(-(x-6+6  )/6A )  ,  (x-6+6.)/6  >  0 

1 

0  ,  elsewhere  , 

which  is  discontinuous  at  the  point  x  ■  6-6^. 


•  *  •*  •'  "*«*  •**  * 

^ >**  .**  s’*  /■  . 

•  v  -.'v-V-'-' 


.s 


>• 


•  ^ •  *  -  *  - s  rc a. iv- ^  ii' k\ ^r^rrrr  < ■  *  ■ 


A  generalization  of  (3.5)  is  straightforward.  let  _t_  “  (tQ*«**»t  )  be  a  non¬ 
decreasing  real  sequence  with  tj^  <  t^+n,  i  ■*  0,...,p-n,  {pg,...,Pn}  be  an  increasing 


integer  sequence,  and  _r_  *  {t.  «  t  },  j  m  0,...,n,  Pg  “  0,  pR  “  p  be  a  subsequence  of  t. 

3  j 


Let  also  at  -  Ti-Tg'  a2  “  T2  "  T1'*",an  “  Tn  “  Tn-1 


and 


b1  '  to* 


'bp  -  tp  -  tp. 


■  t  *■•••  a1  “  b1+...+bp^  *2  “  bpl+l 


+. . .+b 


a  -  b  .+...+b  .  Then  by  the  same  argument  as  in  deriving  (2.5) 

n  p  +i  p 

n- 1  n 


fa1z1,n  +”'+  anzn,n<X_T0> 


=  J,  et'i't<i'1,fb1^--+bi.n+biZ1,n+---+bi+n.1Z„,„<X'T°) 

where  q  *  p  -  n+1,  6n  T  t (i— 1 )  is  the  discrete  B-spline  given  by  (2.4). 


Denote  bv  u  (a  Z.  +...+  a  Z  )  the  1-th  moment  of  the  linear  combination  L 
i  1  1,n  n  n,n 


a  Z,  +...+  a„Z„  _.  1  =  0,1,...  and  assume  the  sequence  x_  is  strictly  increasing.  Then 
i  l ,n  n  u 

the  following  recurrence  formulae  hold. 


W»<aizi  +•••+  a  Z  ) 
1  1  1,n  n  n,n 


(3.6) 


(■:*)•’  ^  -riT)w,.„ .....  v,v,,n' 


where  the  initial  values  are  equal  to 

j+1  j+1 


;j(a1,Z1  ^  -  (x3  -  xjj  )/[( j+1)(x1-xQ)J,  (j  -  0,1,...)  . 


(3.7) 


(3.8) 


I  (?) ' 

T  1  (x  -x  )u  ,(a,Z, 

j-1  3 

1  1  p  m  n-j  1  1,n 

m-0 

,  (p  ■  0,1, « • • ,p  ^  n) 

n 

J-1 

I  (?) 

IT  <vx«,IWai*i,n 

j-i  3 

m-0  r 

_  n  n 

i.Z,  )  -  x  -x 
j  j,n  p  0 


n- 1 


m*p 

,  (p  -  0,1, ...,n-1> 


The  proof  of  formulae  (3.6),  (3.7)  and  (3.8)  follows  applying  Lemma  3.1  to  the 
corresponding  formulae  for  the  moments  of  B-splines  (see  Newman  1980). 
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used. 

To  emphasize  the  sample  size  n  in  the  sequel,  we  shall  also  use  the  notation 
>  Z2  n  >...>  Zn  n  for  the  order  statistics  of  the  sample  Z^ ,Z2, • • • , Zn.  Let  x_  » 
{Xq,x^,.. -»xn}  be  a  non-decreasing  real  sequence  with  x0  <  x^  and  a^  »  x^-xq,  a2  » 
x2-x1,...,an  -  x^x^,. 

LEMMA  3.2. 

Let  n  >  2, 


[X-Xq 

x  -x„  ^a.Z.  .+...+a  ,Z  .  *x 

n  0  1  1,n-1  n-1  n-1,n-1 

x  -x 

+  x  "xn  fai+aozi  ,+..*+a  Z  *x  x0^1 

n  0  12  1,n-1  n  n-1, n-1  I 


V  + 


f*'  >  V'  ^  ^  A^  .v.vv 
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